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Abstract 

Some gamma-ray burst (GRB) spectra exhibit high energy tails 
with the highest photon energy detected at 18 GeV. The spectral slope 
of the high-energy tails is sufficiently flat in vF v to consider the possi- 
bility of their detection at still higher energies. We calculate how many 
bursts can reasonably be expected above a given energy threshold for 
a cosmological distribution of bursts satisfying the observed apparent 
brightness distribution. The crucial point is that the gamma-ray ab- 
sorption by pair production in the intergalactic diffuse radiation field 
eliminates bursts from beyond the gamma-ray horizon t 77 ~ 1, thus 
drastically reducing the number of bursts at high energies. Our results 
are consistent with the non-detection of bursts by current experiments 
in the 100 GeV to 100 TeV energy range. For the earth-bound detec- 
tor array MILAGRO, we predict a maximal GRB rate of ~ 10 events 
per year. The Whipple Observatory can detect, under favorable con- 
ditions, ~ 1 event per year. The event rate for the HEGRA array is 
~ 0.01 per year. Detection of significantly higher rates of bursts would 
severely challenge cosmological burst scenarios. 
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1 Introduction 



The lack of a significant large scale anisotropy in the angular distribution of 
gamma-ray bursts (GRBs) argues in favor of their cosmological origin (e.g., 
Meegan et al. 1995; Briggs et al. 1996; Tegmark et al. 1996). Assuming 
they are all standard candles, the observed number of bursts at a given flux 
relates directly to a number of sources at a given redshift. The maximum 
redshift sampled by BATSE under these assumptions is z max ~ 2 (e.g., Co- 
hen & Piran 1995). Evolution of the burst population can modify this value. 
If bursts would happen to be more active in recent cosmic history, the max- 
imum redshift would be lower. However, the maximum redshift cannot be 
much lower than ~ 1 in any scenario in which the bursts trace the distri- 
bution of galaxies. Galaxy clusters are known to be concentrated toward 
the supergalactic plane (e.g., Tully 1992; Kolatt, Dekel, & Lahav 1995). 
The brightest (and thus nearest) bursts would reflect the mass concentra- 
tions of the nearby local universe. However, no significant deviation from 
isotropy has been found in the 3B-catalog data corrected for sky exposure 
(Hartmann, Briggs, & Mannheim 1996). 

If GRBs are indeed seen to a redshift of ~2, absorption of the high en- 
ergy gamma-rays by diffuse background radiation must become important 
at photon energies greater than 30 GeV. Since gamma-ray emitting blazars 
are also subject to pair absorption, and since their redshifts are often accu- 
rately known, one can determine the density of the diffuse radiation fields 
from the infrared to the UV by measuring their gamma-ray cutoff energies 
(Stecker, de Jager, & Salamon 1992; de Jager, Stecker, & Salamon 1994; 
Dwek & Slavin 1994; Mannheim et al. 1996; Madau & Phinney 1996). For 
blazars, detailed spectral models exist which can discriminate between in- 
ternal and external absorption (e.g., Blandford & Levinson 1995; Mannheim 
1993). 

In fact, one could calibrate the distance scale to gamma-ray bursts by 
measuring their cutoff energies due to external absorption using the back- 
ground radiation density inferred from blazar observations. For such a pro- 
gram to be successful, nature must provide burst spectra reaching very high 
photon energies. Currently, our knowledge about high-energy burst spec- 
tra is rather sparse, and it is not known whether they exhibit an intrinsic 
turnover in the 10GeV~30TeV region or not (Hurley 1996). However, 
burst statistics are consistent with all BATSE bursts having high-energy 
tails such as the ones observed by EGRET (Dingus et al. 1995). Some of 
the bright GRBs which received exposure by EGRET were detected at GeV 
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energies without evidence for a spectral rollover (Hurley et al. 1994). This 
raises the hope that the intrinsic spectra continue to still higher energies and 
that one could in turn use the expected exponential turnover to set a dis- 
tance scale. In fact, the relativistic cosmological fireball model of Meszaros &; 
Rees (1993) predicts that the GRB spectra have high-energy tails in general 
(Meszaros, Rees, & Papathanassiou 1994). Nevertheless, numerous array 
and Cerenkov telescope burst searches at high energies have only provided 
upper limits; e.g. Kieda et al. (1996) (CASA-MIA), Alexandreas et al. 
(1994) (CYGNUS), Aglietta et al. (1992) (EAS-TOP), Krawczynski et al. 
(1996), Funk et al. (1996) (HEGRA), Allen et al. (1995) (MILAGRO), 
Amenomori et al. (1995) (Tibet), and Connaughton et al. (1995) (Whip- 
pie). 

The purpose of this work is to provide quantitative predictions of the 
expected burst rate as a function of photon energy accounting for cosmic 
absorption and to compare the result with experimental limits. We will 
start with the computation of the gamma-ray horizon and show templates 
of absorbed spectra. In a second step, we use the gamma-ray horizon to 
determine the number of bursts at a given detection threshold. Finally, 
we determine the limiting sensitivity required for a detection at a given 
threshold gamma-ray energy and compare this with the sensitivity of current 
experiments. 



2 Gamma-ray absorption from 10 GeV to 100 TeV 

Gamma-rays of energy E propagating from a distant source at redshift z to- 
ward a terrestrial observer can be absorbed by inelastic interactions with low 
energy photons of present-day energy e from an isotropic diffuse background 
radiation field. The dominant process is pair creation + 7e — ► e + + e~ . 
The threshold energy for this process is given by 

2(m c c 2 ) 2 

£th " £7(1-^(1 + 3)2' [L) 

where [i = cos 9 denotes the cosine of the scattering angle. Since the soft 
photon density varies strongly with energy for typical radiation fields, the 
optical depth r 77 must also vary reflecting the number density of target 
photons at the resonant energy oc E^ 1 . The pair creation cross section is 
given by 

^ = ^(1 - P 2 ) [W 2 - 2) + (3 - /3 4 ) In (I±|)" 



cm 2 , (2) 



2 



where {3 = y/1 — I/7 2 with 7 2 = e/e t h, and where <tt denotes the Thomson 
cross section. For the computation of the optical depth we use the geodesic 
radial displacement function dl/dz = ^-[(1 + z)E(z)]~ 1 where E{z) is given 
by equation (13.3) in Peebles (1993). For a cosmological model with = 1 
and A = 0, the function E{z) simplifies to (1 + z) 3 / 2 . Hence, we obtain the 
optical depth 
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for a nonevolving present-day background density nt>, i.e., n' h {z,e')de' = 
(1 + z) 3 rtb(e)<ie, where the dash indicates comoving frame quantities. At 
some time in the past, the background photon density was built up by 
bursts of massive star formation during the era of galaxy formation. As a 
result, the photon density evolution should be more shallow than oc (1 + z) 3 
at this epoch, and the corresponding 7-ray absorption should be somewhat 
weaker. We will discuss this effect only qualitatively below. The shape of 
the present-day diffuse background density nt,(e) is obtained by averaging 
over various galaxy formation models presented by MacMinn & Primack 
(1996). We multiplied this shape by a small factor to obtain agreement 
with the background density (including the contribution from the 2.7 K 
microwave background) estimated by Beichman & Helou (1992) in the far- 
infrared (FIR) (for a modest galaxy luminosity or density evolution 7 = 2), 
and by Madau & Phinney (1996) in the near- infrared (NIR) through UV 
range (Fig.l). The predicted background densities depend sensitively on 
which density or luminosity evolution and maximum redshift are assumed 
and can vary by an order of magnitude. Hence, it is very important to 
tighten existing constraints on the actual background density (Biller et al. 
1995) by obtaining better galaxy luminosity functions from deep galaxy 
surveys or by direct measurements of blazar gamma-ray (external!) turnover 
detections. 

We integrate numerically the optical depth function and solve for the 
gamma-ray horizon r~ n {E , z Q ) = 1. Results are shown in Fig. 2 for two 
values of the Hubble constant. Templates of absorbed infinite power law 
spectra are shown in Fig. 3 for various redshifts. 

Note that the adopted diffuse radiation background does not have a 
strong dust component which would be responsible for TeV absorption of 
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Mrk421, as suggested by de Jager, Stecker, & Salamon (1994). Theoretical 
modeling of the broadband continuum spectrum of Mrk421 over 18 orders 
of magnitude in frequency predicts a cutoff due to internal gamma-ray ab- 
sorption for Mrk421 (Mannheim et al. 1996) which could explain the lack 
of gamma-rays above TeV observed by Whipple. This is in contrast to the 
assumption by de Jager et al. (1994) of a straight power law for the in- 
trinsic spectrum extending well beyond TeV. Further observations will show 
whether the cutoff is variable with flux as expected from the intrinsic model 
or whether it is fixed as expected from absorption by an anomalously high 
external dust component. 



3 Expected total number of bursts 

In our straw person's model, the bursts represent standard candles out to a 
redshift of z max = 2.1 (Cohen & Piran 1995). Assuming no source evolution 
and a galaxy density of ~ 1CT 2 /i 3 Mpc~ 3 (H Q = WOh km s _1 Mpc -1 ), the 
observed BATSE brightness distribution can be matched if ~ 10~ 6 bursts 
(each of energy L a = 7 x 10 50 /i -2 ergs) occur per year per galaxy. The 
number of bursts above a gamma-ray energy E is then given by the number 
of bursts in the cosmic volume enclosed by the gamma-ray horizon at that 
energy, i.e. the hyperplane T^(z , E ) = 1. The absorbed energy from 
bursts beyond the horizon is re-emitted in other wave bands depending on 
the intergalactic magnetic field (Aharonian, Coppi, Sz Vdlk 1994), possibly 
with time delays indicative of intergalactic magnetic fields (Plaga 1995). 
For burst spectra as steep as those observed with EGRET in the GeV range 
(Hurley et al. 1994), the "pileup" of absorbed energy below E a is negligible. 
The number of bursts per year at a given threshold energy E Q and for a 
volume burst rate n Q is then given by 

N{ >E^^(£f j^ 2j 0— (4) 

[see equation (13.61) of Peebles (1993), with an extra factor (1+z) account- 
ing for the redshift of the burst rate] . The integral simplifies considerably 
for = 1 and A = 0. In this case, the angular size distance y(z) is given by 
y(z) = 2[1 - (1 + z)~ 1/2 } and E(z) = (1 + zfl 2 yielding 
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For the volume burst rate n a = 6.9 x 10 9 /i 3 Mpc 3 yr 1 and a BATSE trig- 
ger efficiency of -/Vbatse/-^ = 0.3, we obtain JVbatSE = 300 yr -1 integrating 
to 

z o — z ma,x- The dependence on the Hubble constant enters equation (^) 
only through the value of z = z (E , H Q ) and is therefore rather weak. Nev- 
ertheless, it is an interesting possibility to determine the Hubble constant 
by number counts of GRBs. 

In Fig. 4 we show the number of bursts versus detection threshold demon- 
strating the rapid shrinking of the observable universe with increasing E 
due to pair absorption. At energies of ~10 TeV, the standard candle scenario 
predicts 2 ~ 4 (h = 0.5 ~ 0.75) bursts per year. 



4 Limiting sensitivity 

The number of bursts which can be detected above a given energy threshold 
depends also on the particular instrumental sensitivity F , the spectral shape 
and the duration of the bursts. Suppose we want to measure all GRBs 
above an energy E (within redshifts z < z ), then what is the necessary 
flux sensitivity of our detector to achieve this goal? In the straw person's 
model, this problem is easy to solve. For a given threshold E a the observed 
flux obeys the relation 

L[>E (l + z )} 
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again referring to Peeble's (1993) notation. The standard candle luminosity 
normalized for a BATSE burst duration 5t is given by 

L[> Eo (l + z )] = 7 x 10 5( V 2 ~" (|i) 2 ~ Q (1 + z f~« , (7) 

adopting a high-energy power law spectrum with photon index a > 2 above 
an intrinsic break energy E^ ~ 1 MeV. With = 1 and A = 0, this corre- 
sponds to 



!i + Z 



1— a / X+ \ -1 / TP \ 2-a 



for the source flux at the distance of the gamma ray horizon. All other 
sources will be within the horizon, thus closer, and thus brighter. Although 
the burst luminosity in the standard-candle picture is fixed (5t ~ 1 s), 
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we leave the possibility open that the high-energy tails persist longer than 
the typical BATSE bursts. Hurley et al. (1994) indicate spectral slopes 
a = 2.2 ~ 3.7 and extended or delayed durations of 30s ~ 90 minutes for 
EGRET detected bursts. 

In Fig. 5 we plot the limiting sensitivity for a = 2.2 and a = 2.7 and 
for high-energy tail durations extended relative to the BATSE bursts by 
factors 1 and 100. For comparison, we show current experimental limits in 
the same plot demonstrating that they can in principle detect all bursts. 
The apparently paradoxical increase of the limiting flux with energy reflects 
the rapid shrinking of the gamma-ray horizon: The higher the energy, the 
closer, and thus the brighter, the bursts from the gamma-ray horizon are. 
More pessimistically, if the intrinsic spectra were steeper than the EGRET 
bursts, the sensitivity of current experiments is insufficient to detect all of 
them. In view of the very small total burst rates at high energies (see Fig.4), 
this would practically comply with zero expected event rates. 

5 Conclusions 

We estimated the expected gamma-ray burst rate as a function of threshold 
energy. The straw person's scenario considers standard candles at cosmo- 
logical distances emitting an intrinsically unabsorbed high-energy power law 
spectrum. With these simplifying assumptions, we predict a maximum of 
20 ~ 40 bursts per year in the TeV range, roughly an order of magnitude 
less above 10 TeV, and practically none at 100 TeV. The perhaps more real- 
istic case of a steep power law luminosity distribution of the bursts does not 
significantly affect the maximum redshift z max of their cosmological distri- 
bution (Cohen & Piran 1995), and therefore it does not change the expected 
burst rates. An enhancement of the burst rate in the 10 ~ 100 GeV range 
by factors of ~ 2 is possible if evolution of the background density is impor- 
tant at redshifts ~ 1. The sensitivities of current experiments are sufficient 
to detect bursts if their spectra are not steeper than a ~ 2.7. Taking 
into account sky exposure and triggering efficiency (see Tab. 1), we rule 
out the possibility of cosmological burst detection by CYGNUS, EAS-TOP 
and CASA-MIA (and similar experiments with thresholds at ~ 100 TeV). 
Due to its low threshold, MILAGRO or future low-threshold air Cerenkov 
telescopes could at most detect ~ 10 events per year. The Whipple Ob- 
servatory, the Tibet air shower array, and HEGRA have very low expected 
detection rates (0.01 ~ 1 events per year). If the burst delays at high en- 
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Table 1: Gamma-ray burst rates predicted for various experiments (0 = 1, 
A = 0, h = 0.5). Note that a Hubble constant h = 0.75 would increase all 
rates by a factor ~ 2. 
Experiment Trigger Energy GRB Rate 

Efficiency Threshold (TeV) (year -1 ) 



CASA-MIA 


0.05 


100 


0.0003 


CYGNUS 


0.05 


100 


0.0003 


EAS-TOP 


0.05 


100 


0.0003 


HEGRA-Scintillators 


0.05 


25 


0.02 


CASA-DICE 


0.0004 


20 


0.0002 


HEGRA-AIROBICC 


0.006 


13 


0.01 


MILAGRO 


0.05 


0.2 


10 


Tibet 


0.05 


10 


0.1 


Whipple^ 


0.006 


0.3 


1.0 



a in order to start follow-up observations slewing the Cerenkov telescope can take up 
to lh 



ergies are shorter than the telescope slewing time, the expected burst rate 
for air Gerenkov telescopes would be down by ~ 0.01 owing to their small 
field of view. The major uncertainties of the theoretical expectation values 
lie in the cosmic IR-to-UV photon density (~ 10 in the FIR), the Hubble 
constant (factor < 2) and the effect of forward cascading in the intergalactic 
medium which counteracts the pair absorption depending on the intergalac- 
tic magnetic field strength (Protheroe &; Stanev 1993; Aharonian, Coppi, & 
Volk 1994; Plaga 1995). A low background density, a large Hubble constant 
and weak intergalactic magnetic fields would act to increase the high-energy 
burst detection rate. On the other hand, the high infrared background in- 
ferred by de Jager, Stecker & Salamon (1994) and by Dwek & Slavin (1994) 
assuming TeV absorption of Mrk421 would further reduce the number of 
bursts above ~ 1 TeV by a factor of ~ 10. In any cosmological scenario, 
pair absorption by diffuse background radiation leads to a rapidly decreas- 
ing number of bursts above ~ 10 GeV. By contrast, bursts originating in 
a galactic halo would not be affected by absorption which occurs only on 
cosmic scales. Halo bursts should therefore be much more numerous than 
cosmological bursts at high energies, as has been pointed out by Alexandreas 
et al. (1994) and Hurley (1996). In fact, if one clings to the assumption of 
infinite power law tails, the nondetection of GRBs by current experiments 
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is in agreement with the straw person's scenario of a cosmological origin 
of bursts, but it would require intrinsic spectral cutoffs below ~ 100 GeV 
in halo models. On the other hand, just a few burst detections above ~50 
TeV would render cosmological models highly unlikely. If intrinsic cutoffs 
would be important, it is interesting to note that the predictions for halo 
and cosmological burst models are again quite different. In halo models, 
the brighter bursts should have cutoffs at lower energies due to their higher 
compactness, whereas cosmological models predict the opposite, since the 
bright sources would be near sources with less external absorption. This 
crucial test could be performed by future missions that target the favorable 
energy range 10 ~ 100 GeV. 



6 Acknowledgements 

We thank Matthew Baring, Trevor Weekes and our referee Ocker de Jager 
for their comments. This work was supported in part by NASA under grant 
NAGW 5-1578. D.H. and K.M. are grateful for hospitality and support 
(under NSF grant PHY94-07194) during their stay at the ITP, where this 
collaboration began. B.F. is supported by the BMBF, Germany, under 
contract 05 2WT164. K.M. acknowledges travel support by the Deutsche 
Forschungsgemeinschaft under grant Ma 1545/5-1. 



8 



7 References 



Aglietta, M., et al. 1992, II Nuovo Cimento, 15C, 441 
Aharonian, F.A., Coppi, P.S., & Volk, H.J. 1994, ApJ, 423, L5 
Alexandreas, D.E., et al. 1994, ApJ, 426, LI 

Allen, G. E., et al. 1995, Proc. 24th Int. Cosmic Ray Conf. (Rome), 3, 516 
Amenomori, M. et al. 1995, Proc. 24th Int. Cosmic Ray Conf. (Rome), 2, 
112 

Beichman, C.A., & Helou, G. 1991, ApJ, 370, LI 
Biller, S.D., et al. 1995, ApJ, 445, 227 
Blandford, R.D., & Levinson, A. 1995, ApJ, 441, 79 
Briggs, M. S., et al. 1996, ApJ, 459, 40 
Cohen, E., & Piran, T. 1995, ApJ, 444, L25 

Connaughton, V., et al. 1995, Proc. 24th Int. Cosmic Ray Conf. (Rome), 
2, 96 

de Jager, O.C., Stecker, F.W., & Salamon, M.H. 1994, Nature, 369, 294 
Dingus, B., et al. 1995, ApSS, 231, 187 
Dwek, E., & Slavin, J. 1994, ApJ, 436, 696 

Funk, B., et al. 1996, Proc. 3rd Huntsville Symp. on GRBs, ed. C. 

Kouveliotou, M. Briggs, & G.J. Fishman, (New York: AIP), in press 
Hartmann, D., Briggs, M.S., & Mannheim, K. 1996 ApJ, in press 
Hurley, K., et al. 1994, Nature, 372, 652 
Hurley, K., 1996, Space Sci. Rev., 75, 43 

Kieda, D., et al. 1996, Proc. 3rd Huntsville Symp. on GRBs, ed. C. 

Kouveliotou, M. Briggs, & G.J. Fishman, (New York: AIP), in press 
Krawczynski, H., et al. 1996, Proc. 3rd Huntsville Symp. on GRBs, ed. C. 

Kouveliotou, M. Briggs, &: G.J. Fishman, (New York: AIP), in press 
Kolatt, T., Dekel, A., & Lahav, O. 1995, MNRAS, 275, 797 
MacMinn, D., & Primack, J. 1996, Space Sci. Rev., 75, 413 
Madau, P., & Phinney, E.S. 1996, ApJ, 456, 124 
Mannheim, K. 1993, A&A, 263, 267 

Mannheim, K., Westerhoff, S., Meyer, H., & Fink, H.-H. 1996, A&A, in 
press 

Meegan, C. A., et al. 1995, ApJ, 446, L15 
Meszaros, P., & Rees, M.J. 1993, ApJ, 405, 278 
Meszaros, P., Rees, M.J., & Papathanassiou, H. 1994, ApJ, 432, 181 
Peebles, P.J.E. 1993, Principles of Physical Cosmology, (Princeton: Univ. 
Press) 

Plaga, R. 1995, Nature, 374, 430 



9 



Protheroe, R. J., & Stanev, T. 1993, MNRAS, 264, 191 

Stecker, F.W., de Jager, O.C., & Salamon, M.H. 1992, ApJ, 390, L49 

Tegmark, M., Hartmann, D.H., Briggs, M.S., & Meegan, C.A. 1996, ApJ, 

in press 
Tully, R.B. 1992, ApJ, 388, 9 



10 



Figure Captions 



Fig.l: Spectral energy distribution of the diffuse background radiation. 
Components are the 2.7 K cosmic microwave background (CMB), dust and 
stellar light from galaxies. Solid line: FIR-to-UV background adopted in 
the present work (which closely resembles the spectrum obtained by aver- 
aging over the various cold + hot dark matter models of galaxy formation 
of MacMinn & Primack 1996). Straight line segment: FIR background pho- 
ton density inferred by Dwek & Slavin (1994) assuming TeV absorption for 
Mrk421. Limits: Experimental upper limits on the background IR density 
obtained by Biller et al. (1995). 

Fig. 2 Gamma-ray horizon r 77 = 1 for the diffuse background radiation 
shown in Fig.l. Solid line: h = 0.5. Dashed line: h = 1.0. 

Fig. 3 Templates of absorbed a = 2 spectra at various redshifts at which 
blazars have been detected with EGRET (h = 0.5). 

Fig.4 GRB rate as a function of energy threshold. Solid line: h = 0.5. 
Dashed line: h = 1.0. 

Fig.5 Solid lines: Limiting sensitivity for detection of all GRBs as a function 
of threshold energy assuming a high-energy tail spectral index a (h = 0.5). 
Dashed lines: Same for high-energy tails with extended durations relative 
to the BATSE bursts by a factor of 100. Symbols: current experimental flux 
limits for burst detection (typically assuming St ~ 30 s). 
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